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SUMMARY 


The major objective of the present study was to investigate 
the combustion characteristics of liquid hydrazine strands under 
both steady state and oscillatory conditions* A theoretical model 
of the combustion system was developed and compared with the 
experimental results* 

The theoretical model of the system assumed one-dimensional 
flow and used global kinetic parameters to characterize 
the hydrazine decomposition flame* The properties of the liquid 
phase were assumed constant; however, a variable property gas phase 
solution was obtained through the use of a modified Howarth 
transformation. Conventional low pressure phase equilibrium was 
used to characterize the liquid surface* The equations resulting 
from the theoretical model were solved using a perturbation analysis 
with the amplitude of the imposed pressure oscillations taken as 
the perturbation parameter* Thus, in zero order, the model corresponds 
to steady strand combustion while the first order model takes into 
account the effects of a small amplitude, sinusoidal pressure 
disturbance . 

The steady experimental results were obtained using a steady 
strand burner* With this apparatus, steady strand burning rates 
were obtained for mean pressures ranging from 0*32 up to 42 atm. 

Liquid temperature disturbations and surface temperatures were 
measured in the range of 1*0 to 20.5 atm* 

For the unsteady measurements, an oscillatory version of the 
steady strand apparatus was developed* With this apparatus, the 
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burning rate response of hydrazine to imposed pressure oscillations 
was measured. The apparatus permitted the amplitude of the 
oscillating pressure to be varied independently of the frequency 
at a given mean pressure. Oscillatory burning rates were obtained 
in the mean pressure range of 1 to 10 atm with frequencies varying 
from 0.40 to 5.2 Hz and pressure amplitudes up to 35% of the mean 
pressure. 

The experimental steady strand burning rates were found to 
be dependent on the liquid hydrazine purity. Adding water as an 
impurity to the fuel resulted in a reduction in the burning rate. 
However, for the highest purity fuels tested (98.6 and 99.4%) the 
burning rate of the two fuels was not appreciably different. 

The diameter of the sample tube was found to have an effect 
upon the experimental burning rate. As tube size increased, the 
burning rate decreased due to the reduction in the effects of surface 
tension for the larger tubes c By plotting the strand burning rate 
as a function of the reciprocal of the tube diameter, linear plots 
were obtained allowing the determination of the fundamental strand 
burning rate (burning rate in an infinite diameter tube) by 
extrapolation. 

For subatmospheric pressures the fundamental burning rate 
varied as the square root of pressure; for pressures greater than 
atmospheric, the fundamental burning rate varied linearly with 
pressure. The theory was matched to these results by assuming a first 
order decomposition reaction for subatmospheric pressures and a 
second order reaction for pressures greater than atmospheric. 
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The liquid temperature experimental results gave a good 

prediction of the theoretical model assuming constant liquid 

-3 2 

properties and by using a value of 1«.26 x 10 cm /sec for the 
liquid thermal diffusivity of hydrazine „ The liquid surface 
temperature results were also adequately predicted by the theoretical 
model, justifying the low pressure phase equilibrium assumption. 

The steady strand results indicated that the theoretical 
model was in good agreement with all the available data on hydrazine 
strand combustion. Next, the theoretical model was compared to the 
experimental oscillatory burning rate measurements. It was found 
that the oscillatory measurements of both amplitude and phase angle 
of the liquid surface oscillations were in good agreement with n the 
theoretical model if the decomposition reaction was assumed to be 
second order with an activation energy of about 40 kcal (E=15) for 
mean pressures in the range of 1-10 atm. Unlike the steady state 
where activation energy had little effect upon the predicted burning 
rate, for the oscillatory case activation energy had a pronounced 
effect upon both the amplitude and phase angle of the unsteady burning 
rate . 

The theoretical results showed that for activation energies 
characteristic of hydrazine (E=10-15) an increase in the response 
of the combustion process occurred due to interaction with transient 
liquid phase effects. This yielded a band of frequencies where the 
combustion process exerted sufficient amplifying power to provide a 
mechanism for driving combustion instability. As pressure increased 
the amplifying power of the combustion process increased and this 
frequency band moved to higher frequencies. 



CHAPTER I 


INTRODUCTION 

1. 1 General Statement of the Problem 

The problem of combustion instability has plagued liquid 
rocket engine designers for a number of years . Since this phenomenon 
severely reduces the efficiency and thrust of the engine and in 
extreme cases may even cause destruction of the rocket, there is 
much interest in gaining an understanding of the process. Briefly, 
combustion instability results from the coupling of the combustion 
processes with the fluid dynamic processes in a combustion chamber. 

As a result of this coupling, the combustion supplies oscillatory 
energy to sustain any oscillations of the fluid dynamic variables 
(pressure, velocity, etc.). If sufficient damping exists in the 
system, the energy will be dissipated more rapidly than it is supplied 
and the oscillations will decay. Thus, combustion instability can 
be prevented by either increasing the damping of the system or by 
decreasing the coupling of the combustion and fluid dynamic processes. 

The main method of suppressing combustion instability has 
been to use acoustic liners, Helmholtz reasonators and other damping 
devices in the combustion chamber. However, even with these damping 
devices one can never be sure a combustion chamber will be stable 
until test firings are made. The most desirable approach to the 
problem is to gain an understanding of the coupling between the 
combustion and fluid dynamics of the system and then to develop 
means of reducing this effect. 
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There have been basically two approaches to the theoretical 
description of the processes taking place in a combustion chamber. 
One method is the heuristic approach of Crocco [1, 2, 3] or the so- 
called time-lag theory. Crocco postulates the existence of a time- 
lag between the time when an element of fuel is injected into the 
combustion chamber and the time when it is gasified by combustion to 
final products. However, the time-lag can only be deduced from 
experimental rocket firings and no detailed description of the 
coupling between combustion and the gas dynamics is offered. 

The other approach has been due to Priera [4, 5, 6] who 
numerically solves the gas dynamic equations for the combustion 
chamber. Two of the most important parameters of the analysis are 
the thermal energy and gas release rates of the combustion process. 
With Priem*s approach, attention is focussed on the nonsteady 
burning rate response of an individual droplet to a fluctuation in 
the surrounding gaseous flow field (either a pressure or a velocity 
fluctuation or both) • This burning rate response is then used as the 
source term in the gas dynamic equations. 

The above discussion has indicated the importance of the 
burning rate response of a liquid to fluctuations in its environment. 
Thus the subject of the present investigation is the study of the 
burning rate response of a liquid to fluctuations in ambient 
conditions. For ease of interpretation of the results, only 
pressure fluctuations are considered. In addition, the pressure 
fluctuations are assumed small when compared to the mean pressure. 

The particular fuel chosen for study is hydrazine since the hydrazine 
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family of fuels comprise some of tihe most important fuels used in 
current space rocket applications. 

In the next section a review is made of the available 
literature on steady hydrazine combustion. Emphasis is placed on 
hydrazine strand burning since this experimental technique is used 
in the present study. In the following section, the literature on 
liquid combustion instability, in particular the unsteady burning 
rate response of a liquid to fluctuations in its environment, is 
reviewed. 


1.2 . Previous Related Studies 

1.2.1 Steady Hydrazine Combustion 

iNlumerous studies have been conducted on steady state liquid 
hydrazine combustion as droplets and also as liquid strands. The 
literature on hydrazine droplet combustion has been reviewed in 
detail in Reference [7]. By considering the bulk of the data on 
hydrazine droplet combustion. References [7-11], Reference [7] 
concludes that a monopropellant gas phase model is more realistic 
than a bipropellant model for hydrazine. 

The work on hydrazine strand combustion has been mainly 
concerned with the effect of pressure on hydrazine burning rates. 

With the strand burner technique, liquid hydrazine is placed in an 
open-ended test tube which is contained within a windowed pressure 
vessel. After ignition, the liquid surface moves down the tube as 
the fuel is consumed by the combustion process. The rate of propagation 
of the liquid surface (the strand burning rate) is then measured as 


a function of pressure. 
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Adams and Stocks [12] investigated the burning rate of 
liquid hydrazine strands in a nitrogen pressurized vessel over the 
pressure range of 1-45 atm. Tube diameters of 3 and 5 mm were used. 
Burning rates were measured by timing the movement of the liquid 
surface, between marks placed 3 cm apart on the tube surface, with 
a stopwatch. The results of Adams and Stocks [12] show considerable 
scatter. However, they conclude that for pressures up to about 10 
atm the burning rate of liquid hydrazine is proportional to the 
square root of pressure. For higher pressures their data Indicates 
that the burning rate is independent of pressure. They also found 
that above a certain pressure hydrazine failed to burn as a liquid 
strand. This upper extinction pressure was found to be a function 
of the tube diameter and also the concentration of the liquid 
hydrazine. The upper extinction pressure increased with increasing 
tube diameter and increasing liquid fuel concentration. In addition, 
the upper extinction pressure increased by substituting helium for 
the nitrogen pressurizing gas. 

Gray and Kay [13] and Gray, et al., [14] extended the work 
of Adams and Stocks [12] to subatmospheric pressures. They found 
that a lower extinction pressure, below which the liquid failed 
to propagate down the tube, also exists for hydrazine. 

The burning rate of liquid hydrazine was found to be 
proportional to the square root of pressure for pressures as low as 
0.4 atm [13 and 14] in agreement with the results of Adams and 
Stocks [12]. The lower extinction pressure increased with increasing 
tube diameter and also increased with decreasing liquid purity. As 
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the lower extinction pressure was approached a slight bubbling 
of the liquid near the liquid surface was observed. 

Antoine [15] studied the burning rate of liquid hydrazine 
strands in the pressure range of 1-19 atm. The tube diameter and 
liquid purity were found to have a significant effect on the burning 
rate. As tube diameter decreased and/or liquid purity decreased, the 
burning rate decreased significantly at a given pressure. For the 
100% concentration liquid hydrazine burning in the largest diameter 
tube tested (12.7 mm), the burning rate varied linearly with pressure 
over the pressure range tested. The extinction pressures for hydrazine 
were not investigated by Antoine, 

Thus, the previous studies on liquid hydrazine strand 
combustion all Indicate that liquid phase purity and tube diameter 
have a significant effect on the burning rate. The results of Adams 
and Stocks [12], Gray and Kay [13] and Gray, et al. , [14] indicate 
that the burning rate of hydrazine is proportional to the square root 
of pressure in the pressure range of 0.4-10.0 atm whereas the results 
of Antoine [15] indicate that for 100% N 2 H^ the burning rate is 
proportional to pressure in the pressure range of 1-10 atm for the 
12.7 mm tube, 

1.2.2 Combustion Instability 

The literature on liquid propellant combustion instability 
is vast (for example. Reference [16] cites 778 references). Therefore, 
only the most important studies on liquid burning rate response to 
imposed ambient oscillations will be considered here. 
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Experimental studies on nonsteady liquid combustion are 
virtually non-existent. The approach used in the past has been to 
infer the burning rate response from the pressure-time curves from 
actual rocket engine firings. Of course, this approach depends on 
some theoretical model of the combustion-gas dynamic problem, and direct 
experimental verification of any unsteady combustion model is 
impossible with this technique. 

Numerous unsteady liquid combustion theoretical models have 
been proposed. Two of the earliest models are quite similar. 

Williams [17] and Strahle [18, 19] considered longitudinal standing 
wave fluctuations acting on the leading edge of a burning droplet* 

Both assumed the classic Burke-Schumann thin diffusion flame existed 
in the interior of the leading edge boundary layer. Williams 
considered only a flat plate whereas Strahle studied both the stagnation 
point and flat plate flows. 

Both of these studies found only rather flat response curves 
as a function of frequency and, moreover, the peak that does exist 
occurs at such a large frequency (1000-10,000 H^) that no influence 
on instability is expected. However, the fallacy with these treatments 
can be traced to the assumption of constant liquid temperature with 
time. It can be shown that the characteristic thermal wave time in 
the liquid phase is at least an order of magnitude greater than the 
characteristic thermal wave time in the gas phase. [20] Thus, on 
this basis, the gas phase is more logically considered quasi-steady 
compared to the liquid phase, especially at low frequencies. In 
other words, the assumption of constant liquid temperature cannot be 
true even at low frequencies. 
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Heidmaun and Wieber [21] assumed that the gas phase burning 
rate of the fuel was equal to the liquid fuel vaporization rate at 
each instant of time. For low frequencies this is a valid assumption. 
However, as in the investigations of Williams [17] and Strahle 
[18, 19], the thermal wave time effects are neglected and the droplet 
temperature, at any instant of time, is considered uniform but 
varying with time. 

T^en and Sirignano [22] recognized the importance of the 
liquid thermal wave time. They considered a fuel evaporating from 
a flat plate with a reacting gaseous boundary layer and a longitudinally 
oscillating external flow. The gas phase was considered quasi- steady 
compared with the liquid phase which is correct in the limit of 
low frequency oscillations. They found peaks in the response 
function which are directly related to the effects of thermal lag 
in the liquid phase. However, it is questionable that this analysis 
indicates the true combustion gas dynamic coupling mechanism since 
the peaks are not of sufficient magnitude to provide instability in 
an actual rocket. 

1.3 >. Specific Statement of the Problem 

As indicated by the preceding discussion, although a number 
of theories of the unsteady combustion of liquid fuels have been 
developed, the direct experimental verification of any of them has 
not been accomplished up to the present time. In addition, all of 
the theoretical models consider a bipropellant gas phase combustion 
process. Some of the most important fuels used in current space 
rocket applications are liquid hydrazine and its derivatives. Earlier 
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work on the combustion of hydrazine indicates that a monopropellant 
gas phase model is more realistic for the hydrazine fuels. [7] 

However, discrepancies exist concerning the reported dependence of 
the steady burning rate of hydrazine on pressure. [12-15] 

Thus, the specific objectives of the present study were: 

1. Study the steady strand combustion characteristics of 

hydrazine both theoretically and experimentally. 

a. Experimentally determine the steady strand burning 
rate as a function of pressure, liquid concentration 
and tube diameter. 

b. Experimentally measure the liquid temperature 
distribution and surface temperature as a function 
of pressure, 

c. Develop a theoretical model of the system and 
compare directly with the experimental results. 

2. Study the oscillatory combustion characteristics of 

hydrazine both theoretically and experimentally. 

a. Experimentally measure the oscillatory burning rate 
(both magnitude and phase angle with respect to 

to the imposed pressure oscillations) of hydrazine 
as a function of mean pressure and the frequency 
and amplitude of the imposed pressure oscillations. 

b. Develop a theoretical model of the oscillatory 
system and directly compare the theoretical 
predictions with the experimental results. 



CHAPTER IX 


EXPERIMENTAL APPARATUS 

As indicated in the Introduction^ a major objective of the 
present study was to experimentally measure unsteady hydrazine 
burning rates » However, before studying the unsteady combustion 
characteristics of hydrazine, the steady state characteristics were 
studied thoroughly* 

In this chapter the apparatus used to obtain the experimental 
data are described* First, the steady strand combustion apparatus 
is discussed followed by a description of the oscillatory combustion 
apparatus. Since the basic configurations of the two apparatus were 
the same, many of the components used in the apparatus, such as the 
ignitor, pressure gages and instrumentation, were similar. These 
components are described in detail in the section on the steady strand 
apparatus ♦ 

2.1 Steady Strand Apparatus 

The liquid strand combustion apparatus was used to study 
the steady state combustion characteristics of hydrazine* With this 
technique the burning rate of the fuel is determined by measuring 
the constant rate of regression of the liquid surface as the fuel 
is consumed by the combustion process. In addition to measuring 
burning rates, liquid temperatures were measured using a thermocouple* 
In the next section the steady strand burner is described in detail 
followed by a description of the test procedure. 
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2,1,1 General Description 

A sketch of the steady strand apparatus is shown in Figure 1, 
The pressure vessel used for the tests had an inside diameter of 
6,4 cm with an inside length of 28 cm and was rated to a pressure of 
400 atm. The windows of the vessel provided a 2.5 cm diameter 
viewing space. The ignitor coil used to ignite the fuel sample was 
made of spiral wound 24 gage nichrome wire about 30 cm in length. 

The voltage applied to the wire varied from between 55 to 65 volts, 

A Heise pressure gage with a range of 0-750 psia was used to measure 
the total pressure in the test chamber. 

Three sizes of pyrex sample tubes were used — 4* 8 and 12 mm 
I,D, each with a wall thickness of 1 mm. The thermocouple was 
located 4 cm from the lower edge of the ignitor coil and 3 cm from 
the bottom of the tube to insure that the combustion wave was steady 
when the temperature measurements were made. 

The optical system consisted of a background light and high 
speed motion picture camera. A PEK Labs. Model No. 911 mercury arc 
lamp was used to provide the back light to the liquid column. The 
shadowgraphs were taken with a Photosonics Model No, 16 mm- lb 16 mm 
camera using Kodak Plus X Reversal film operating at speeds up to 
100 frames per second. The distance measurements from the film 
records were obtained using a Vanguard Motion Analyzer. The Vanguard 
provided a gain of about 25:1 which was adequate for accurate 
resolution of the film records. 

The thermocouples used to record the liquid temperatures were 
constructed with 0.0003 inch 0,D, platinum-platinum 10% rhodium butt 
welded wire. The procedure used to construct the thermocouples is 




Figure 1 Sketch of the Steady Strand Burner Apparatus 
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described in Reference [23]. The thermocouples were stretched 
horizontally (to minimize conduction errors) through holes burned 
in the glass tube and sealed in place with epoxy. In order to 
maintain the reference junction of the thermocouple at a fixed 
temperature, the thermocouple leads were carried outside of the 
test chamber such that the reference temperature was maintained at 
room temperature. 

The thermocouple output was recorded on a CEC Type 5-124 
oscillograph having a flat frequency response to 2000 Hz. The film 
and temperature records were synchronized by a switch closure which 
deflected a galvanometer on the oscillograph and started a light 
streak on the film. A 100 Ez timing signal was placed on the 
oscillograph using its internal timing generator and also on the film 
record using a 100 Hz Wollensak pulse generator. 

Hydrazine (95+% purity) obtained from the Eastman Organic 
Chemicals Company was used in the bulk of the testing. A gas 
chromatographic analysis performed on a sample of this fuel by Bell 
Aerospace Company indicated that it was composed of 98.5% hydrazine, 
1*3% water and Od trace impurities* Some limited testing was also 
done using hydrazine of 99.4% obtained from the Matheson, Coleman 
and Bell Corporation. 

2*1.2 Test Procedure 

The general experimental procedure was similar to that 
employed in References [24-26]* The liquid fuel was placed in a 
glass tube contained within a windowed chamber. After pressurizing 
the chamber with nitrogen to the desired test pressure, the fuel 
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was ignited with a heater coil. Following ignition, the liquid 
burned down the tube with a constant rate, after a short development 
period for the combustion wave. As the liquid surface propagated 
past the window in the chamber^ the rate of regression of the liquid 
column as well as the position of the thermocouple in the liquid 
phase was determined from motion picture shadowgraphs taken through 
the windows of the chamber. Thus the test data consisted of a complete 
liquid phase temperature record as well as the burning rate of the 
fuel. 

Prior to testing on any given day with the steady strand 
combustion apparatus, the thermocouple and optical systems were 
calibrated. The oscillograph was calibrated by applying a known 
voltage to the input of the galvanometer and noting the output 
deflection of the light beam on the trace. The optical system was 
calibrated by positioning a wire of known diameter (measured with a 
micrometer) at the thermocouple location in the viewing space of the 
camera and photographing it. 

After calibration, the following steps were followed during 
a typical experimental run. First the tube and test chamber 
windows were cleaned. The tube was then filled with fuel and 
positioned in the test chamber. The camera was focussed and the 
chamber pressurized with nitrogen to the desired test pressure. The 
ignitor was turned on and then off following ignition; the oscillograph 
and camera were turned on. As the liquid surface came into the 
viewing space of the camera, the synchronization switch was activated 
to synchronize the film and temperature records. Following 



14 


completion of the test all systems were turned off, and the chamber 
was flushed with nitrogen to remove the combustion product gases. 

2.2 Oscillatory Strand Combustion Apparatus 

In order to study the unsteady combustion of hydrazine an 
oscillatory version of the steady strand apparatus was developed. 

As noted previously, in the steady strand combustion case the liquid 
surface regresses at a constant rate down the tube. However, under 
unsteady conditions the surface regression rate is no longer steady and, 
in particular, under oscillatory conditions the surface regression 
rate oscillates about some mean burning rate. The next four sections 
give a detailed description of the oscillatory apparatus and its 
operation. 

2.2.1 General Description 

The experimental apparatus was required to have the 
capability of varying the mean pressure and the amplitude and 
frequency of the pressure oscillations and measuring the response 
of a burning liquid to these parameters. Figure 2 is a sketch of the 
oscillatory combustion apparatus with the above capabilities. 

The major components of this apparatus are the glass tube filled with 
liquid fuel, the camera- lens optical system used to measure the 
fluctuations of the burning liquid surface in response to the imposed 
pressure oscillations, and the rotary valve arrangement used to 
provide the oscillating pressure variations. 

In order to minimize the total flow rate of air required to 
produce a given pressure oscillation, the chamber volume was kept to a 



Figure 2 Sketch of the Oscillatory Strand Combustion Apparatus 
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minimum. The test chamber used in the oscillatory experiments had 
an inside diameter of 2.6 cm with an inside length of 17 cm and 
was leak checked to a pressure of 15 atm. The windows in the 
chamber had a 2.5 cm viewing space. 

A photograph of the test chamber is shown in Figure 3. 

Also shown in the photograph are the camera and lens used to 
photograph the liquid surface. These items are discussed in detail 
in the next section. 

The ignitor coil used to ignite the fuel sample was of 
similar construction to that used in the steady strand burner. The 
static pressure gages sketched in Figure 2 were 0-300 psia Heise 
gages. 

In contrast to the steady strand experiments, only one 
size sample tube was used in the oscillatory tests. The tube size 
used was 8 mm I.D. with a tube wall thickness of about 1 mm. As 
in the steady strand work, the oscillatory burning rate measurements 
were made at a position about A cm from the lower edge of the ignitor 
coil and about 3 cm from the bottom edge to insure that end effects 
were minimal. The 98.6% purity hydrazine obtained from Eastman 
Organic Chemicals Company was used throughout the oscillatory 
testing. 

2.2.2 Photographic System 

Quasi-steady analysis indicated that typical liquid surface 
oscillation amplitudes were on the order of 0.01 mm. In order to 
optically resolve such small distances, the camera-lens system 
sketched on Figure 2 and shown in the photograph in Figure 3 was used 
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Figure 3 Photograph of the Oscillatory Strand Test Chamber 
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in conjunction with a Vanguard motion analyzer. The camera was a 
35 nun Dumont Type 2582 strip camera and the lens had a focal length 
of 150 mm and was 6 mm in diameter. The camera was modified by 
removing its lens and replacing it with an 0.02 inch slot. 

Both the lens and camera were mounted on movable bases so 
that the overall gain of the camera- lens system could be caried. 
Primary gains as high as 5:1 were obtained with this system. 

This gain coupled with the gain of about 25:1 of the Vanguard Motion 
Analyzer provided overall gains of 125:1 which was adequate for 
accurate resolution of the film records. 

Kodak Plus^X Pan film was used throughout the testing. The 
camera speed was varied depending on the test conditions from about 
133 to 400 inches per minute. Timing marks were placed on the films 
using an Adtrol timing pulse generator. 

The background light was supplied by a mercury arc lamp 
identical to the one used in the steady strand work. The background 
light was dispersed by passing it through a ground glass filter. The 
intensity of the light was reduced by removing the focussing mirror 
from the lamp and also by passing the light through a blue filter, 

2*2.3 Pressure Supply System 

The oscillating pressure in the test chamber was established 
using the rotary valve arrangement shown in Figure 2, This method 
is similar to the system described in Reference [27] for establishing 
an oscillating propane gas flame. With this technique an oscillating 
pressure was set up in the test chamber by an oscillating air stream, 
A stream of air was passed through a needle valve and ball valve 
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mounted in parallel. The amplitude of the oscillatory pressure was 
varied independently of the frequency by varying the relative amounts 
of air passing through the two valves. The frequency of the 
oscillating pressure was varied by changing the speed of the DC 
motor used to rotate the ball valve. The mean pressure in the 
chamber was varied by the setting of the needle valve downstream from 
the ball valve o A photograph of the rotary valve arrangement and 
DC motor are shown in Figure The speed of the DC motor was 
reduced about 6:1 with the pulley-V belt system shown in Figure 4. 

Two configurations for the rotary valve system were tested. 
The first configuration tested had the rotary valve system upstream 
of the test chamber. However, by placing the rotary valve system 
downstream of the test chamber, as sketched in Figure 2, a more 
nearly sinusoidal pressure trace was obtained. 

The pressure in the chamber was recorded using a Kistler 
Type 603A pressure transducer with the output displayed on both an 
oscilloscope and an oscillograph for retention. The output from 
the transducer was fed through a Kistler Model 504 Electrostatic 
Charge Amplifier and then to the oscilloscope. From the charge 
amplifier the signal was fed through a CEC Type 1-165 DC amplifier 
to the CEC Type 5-124 oscillograph. The gain of the pressure 
recording system with the oscillograph varied from about 0.8 psi/ 
inch up to about 3 psi/inch. 

The incoming air stream was supplied by a 350 psi shop 
oil-free compressor. Three-eighths inch stainless steel tubing was 
used for all the piping in the pressure supply system. 
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Figure 4 Photograph of the Rotary Valve Arrangement 



Any phase shift between the imposed pressure oscillations 


and the resulting burning rate oscillations were determined from the 
films and the pressure trace registered on the oscillograph. The 
film and pressure records were synchronized by a switch closure 
which deflected a galvanometer on the pressure record and started a 
light streak on the film. Since the camera had no internal lamp, the 
light streak was supplied by an external neon lamp* Fiber optics 
were used to pipe the light from the neon lamp through the camera 
housing to the film. Timing marks were placed on the pressure 
record using the internal timing generator of the oscillograph 
and on the film records by an Adrol timing generator. 

2.2.4 Test Procedure 

The Kistler pressure transducer was calibrated using a dead 
weight tester and the calibration was checked each day of testing with 
a mercury manometer. The optical system was calibrated at the 
beginning of each day of testing and also after each adjustment of 
the gain by photographing the resulting movement of the liquid 
surface when adding a knoxm volume of liquid from a Hamilton Type 
1001-N micro-syringe. For the optical calibration, air was passed 
through the system at about the total flow rate used in testing so 
any effects of the air flow could be taken into account. 

After the daily calibrations, the windows and lens were 
throoughly cleaned. The sample tube was filled with fuel and 
positioned in the test chamber. The test condition in the chamber 
was set by the pressure supply system. The fuel was then ignited 
with the heater coil. Both the camera and oscillograph were turned 
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on; the switch was activated which started the light streak on the 
film and deflected the galvanometer on the pressure trace. After 
completion of the test, the camera and oscillograph were turned off. 
After flushing the combustion product gases from the chamber, the 
tube was removed and cleaned, and the windows and lens were cleaned 
in preparation for the next test. 



CHAPTER III 


THEORETICAL CONSIDERATIONS 

As indicated in the Introduction, this study was concerned 
with the influence of pressure transients on the combustion of a liquid 
monopropellant fuel^ The combustion system which the theoretical model 
must describe consists of a fuel evaporating from a liquid surface and 
then subsequently undergoing exothermic chemical decomposition in the 
gas phase. A portion of the heat evolved in the decomposition 
process flows back to the liquid surface providing the energy 
for further evaporation of the liquid; the remainder of the energy 
released is convected downstream by the gas flow. The main thrust 
of this study was to gain an understanding of the coupling between 
an oscillatory pressure wave and the fuel burning rate. 

In this chapter the theoretical model of the raonopropellant 
combustion process is described. After reviewing some related solid 
propellant theories, the model and assumptions are discussed, 
followed by the development of the equations describing the system. 
These equations are placed in dimensionless form and a perturbation 
solution sought with the amplitude of the oscillating pressure used 
as the small perturbation parameter. 

3.1 Existing Theories 

Since the process is assumed to be one dimensional, some 
insight into the modelling of the experiment can be gained by studying 
solid propellant combustion instability theories. In particular. 
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neglecting condensed phase reactions, solid propellant models are 
applicable to the present problem with modified boundary conditions 
at the two phase interface and slight modifications of the gas phase 
analysis. However, even with the one dimensional assumption the 
equations describing the complete unsteady process are extremely 
complex since they must account for the interaction between the time 
dependent chemical kinetics and gas dynamics. Although the present 
theoretical model is of the one-dimensional experiment, the results 
should lead to a better understanding of liquid response and thus to 
improved droplet response models. 

In contrast to liquid propellant instability studies, solid 
propellant investigators have recognized the importance of the 
condensed phase characteristic thermal wave time for many years. Hart 
and KcCiure [28] discuss the relevant characteristic times for the 
various time dependent processes associated with unsteady solid 
propellant combustion. Based on this study, the various theoretical 
models [29-31] have all assumed a quasi-steady gas phase and completely 
unsteady condensed phase. An exception to this is the recent study 
of T'ien [32] who relaxed the quasi-steady gas phase assumption and 
numerically obtained solutions for an unsteady gas phase as well as 
an unsteady liquid phase. 

However, even with the quasi-steady gas phase assumption, 
analytical solutions cannot be obtained unless further simplifying 
assumptions are made. The need for further simplification arises 
because of the highly nonlinear chemical kinetic terms appearing in 
the quasi-steady gas phase equations. The form of these simplifying 
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assumptions is where the various theories differ* For example^ Krier> 
et al-, [29] assume a zero activation energy for the global form of 
the kinetic term whereas Friedly and Petersen [30] effectively limit 
their solution to infinite activation energies (infinitely thin 
flames) by using the ZelMovich [33] flame theory, Denison and 
Baum [31] use the results of adiabatic laminar flame theory to 
formulate the nonadiabatic solid combustion problem in analytically 
tractable form. Reference [34] offers a more complete discussion of 
the various solid combustion theories. 

The most rigorous approach to the solid propellant combustion 
problem is that due to xHen [32] who models the gaseous reaction 
with global kinetic parameters and numerically integrates the steady 
state gas phase equations. An approach similar to that of T’ien was 
used in the present investigation. The major differences between the 
present analysis and that of T'ien involve phase equilibrium through 
the Clausius-Clapeyron equation and variable gas phase properties 
in the present case as opposed to an Arhenius surface gasification 
rate and constant gas phase properties for T^ien's analysis. In 
addition, T^ien numerically integrated his equations from the surface 
to some fictitious flame length. The present solution method avoids 
the introduction of a flame length by using asymptotic analysis to 
generate starting values for large distances from the surface and 
numerically integrating to the surface, 

3.2 General Model and Assumptions 

A sketch of the theoretical model is shovm on Figure 5. 

The coordinate system chosen is inertial with respect to the origin 
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which is located at the mean position of the burning surface. [35, 36] 

Assuming imposed sinusoidal pressure oscillations, the instantaneous 

position of the liquid surface, x , varies as shown in Figure 5 to 

s 

first order in pressure amplitude. The general appearance of the 
steady state temperature, fuel mass fraction and gas phase reaction 
rate profiles are also sketched on Figure 5. 

For generality, the gas phase transient effects have been 
included in the analysis, similar to the approach used by T*ien 
for solid propellants. [32] The effects of variable properties are 
also included through the introduction of a modified Howarth 
transformation . 

The major assumptions used in the analysis are as follows: 

1. The flow is one-dimensional with a Mach number much 
less than unity and all body forces are negligible. 

2. The flame is premixed and laminar. A one-step, 
irreversible chemical reaction takes place in the gas 
phase and any time lags associated with the chemical 
reaction itself are negligible. 

3. Radiation heat transfer is negligible. An estimation 

of the radiation contribution indicates that it accounts 
for less than 1% of the total energy required for the 
vaporization of hydrazine under steady state 
conditions (see Appendix A) . This is primarily due to 
the fact that hydrazine has a low adiabatic flame 
temperature. 

4. With regard to thermodynamic properties, all gas 
phase diffusion coefficients are equal, all molecular 
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weights are equal and constant, all gas phase specific 
heats are equal and constant, the gas phase mixture 
thermal conductivity varies linearly with temperature 
but is independent of composition, and the liquid is 
composed of a single chemical substance having constant 
properties. 

5. Since the pressure levels are low, the combustion 
products are taken to be insoluble in the liquid 
phase. Chemical reaction is also neglected in the 
liquid phase. 

6. Based on the findings of the steady strand portion of 
the study, the Lewis number is assumed to be unity. 

7. The gas phase is taken to be an ideal gas and at the 
liquid surface fuel mass fraction is related to the 
surface temperature through the Clausius-Clapeyron 
equation. The justification of the use of the 
Clausius-Clapeyron equation is presented in Appendix A. 

8. Since the Mach number is much less than unity, the 
inertial and viscous terms in the momentum equation 
are neglected. 

9. The complete unsteady gas phase equations are used 

but the wavelength of any periodic pressure disturbance 
is assumed to be long compared to the distance from 
the outer edge of the reaction zone to the liquid 
surface. Estimation indicates that this assumption 
is valid for frequencies up to about 50,000 Hz (see 
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Appendix A) . Assumptions 8 and 9 indicate that 
pressure is only a function of time. 

Thus, the theoretical model includes the effects of time 
dependent heat and mass transfer in both the liquid and gas phases. 
Variable property effects are also included in the gas phase and 
global kinetic parameters are used to characterize the hydrazine 
decomposition flame. 

3, 3 Governing Equations 

3.3.1 Dimensional Equations 

With the above assumptions, the equations of overall 
continuity, species conservation and energy in the gas phase and 
overall continuity and energy in the liquid phase were obtained from 
the general equations given by Williams [37]. 

The equations are as follows: 

Gas Phase, x *(t*) < 

7 S 

Conservation of Mass 




( 1 ) 


Conservation of Species 
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Since the gas is ideal 


used; 


= p*RT*/M (4) 

In Equations (2 and 3) the following definitions have been 



(5) 

( 6 ) 

(7) 


Here w* is the rate of reaction of the fuel with a pre-exponential 

iJ 

factor of B*T* and an activation energy of is the heat of 

combustion and is related to the stoichiometric coefficients. 

1 


Liquid Phase, < x *(t*) 

s 

Conservation of Mass 


v| = constant 


Conservation of Energy 


If 

dX^ 


= 0 


( 8 ) 


(9) 


Since only oscillatory solutions are required, the initial 
conditions are irrelevant. The boundary conditions applicable to 
the present problem can be summarized as follows. At the cold end of 
the liquid propellant, the temperature must be constant with respect to 


both X* and t*. 
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X * -.00 

constant) (10) 

At the liquid surface, x*(t*), six conditions must be satisfied* 

s 

The first two are obvious and are mentioned only for completeness. 
That is, the temperature and pressure must be continuous across the 
surface. 


= T* 

s- s+ 


s- s+ 


( 11 ) 

( 12 ) 


Conservation of mass applied across the moving surface reduces to: 


P^(v* - X*) = p*(v* - X*) (13) 


and conservation of energy at the surface yields: 


J s— L J s+ 


A* 

I 


Since the products of combustion are assumed to be insoluble in the 
liquid phase, the gradient of the fuel mass fraction is related to 
the fuel mass fraction itself at the liquid surface as follows: 

r 3Y, 


p*D 


3x* 


s+ 




(15) 


The Clausius-Clapeyron equation relates the fuel mass fraction at the 
surface to the surface temperature 


a* 

'f.s = 


exp 


\-ji] 


(16) 
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where L* is the heat of vaporization of the fuel and a* is a constant 
with the dimensions of pressure. The remaining boundary conditions 
relate the variation of fuel mass fraction and temperature far from 
the liquid surface. Far from the surface, the gas phase chemical 
reaction must go to completion since eventually all the fuel will 
react . Thus 


“>■ 00 

Yp ^ 0 (17) 

for all time. In addition, the temperature must become independent 
of X* since the only energy source, the chemical reaction, has gone 
to completion and heat conduction will smooth out any temperature 
variations with respect to distance. However, the temperature will 
still vary with time since the pressure is a function of time. 
Therefore: 


X* 00 


T*(x*, t^) T*(t^) (18) 


where is a known function. 


3.3.2 Nondimensxonal Equations 

Equations (1-18) are placed in dimensionless form by 
introducing the following variables. 


P = 


P* 

*^00Q 


T = 


V = 


p = 


p* 

p* 

o 


A - 


A^ 

A* 

OOQ 


t = 


t*p* 

OOq 


2 

C V* 
p 


X* 

OOQ 



33 


X = 


x*p* C V* 
p 

X* 


(19) 


The subscript “o signifies a zero order (steady state) quantity 
evaluated at infinity. 

Treatment of variable gas phase properties is accomplished 
by the Howarth Transformation appropriate to the present coordinate 
system by defining 


C V* 

COq 
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p*dx* 


( 20 ) 


x*(t*) 


and by defining a modified mass flux as 


X 


r = pv + ^ / pdx 


I 


( 21 ) 


Xs(t) 

The Shvab-Zeldovich variable 6, is also introduced as a 
new dependent variable. 


e = qYp + T 


( 22 ) 


Substituting Equations (19-22) into Equations (1-4) then 

yields; 


Gas Phase, 0 < r| < " 




(23) 

(24) 

(25) 


P = pT 
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In deriving Equations (24 and 25) use has been made of the 

fact that = -I. The dimensionless reaction rate is defined as 
F 


w 


, A_ ^1+6-n pn-1 


(27) 


where E is the dimensionless activation energy, 


E = 


E* 

RT* 

OOQ 


q is the dimensionless heat of combustion 


- q* 
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(29) 


p ®o 

and A is related to the pre”exponential factor and zero order mass 
flux 
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Also, noting that 


X* T* 


X* 

COQ 


<x>o 


(31) 


and through the use of the ideal gas equation of state, it follows 
that 


pX « P 


( 32 ) 


The dimensionless liquid phase continuity and energy 
equations are 


Liquid Phase, 


< n < 0 
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r = r(t) 


3T ^ 3T 


8^T 


= 0 


(33) 

(34) 


The dimensionless parameter is a constant defined as 

P 0 ^ 0 

6 = P (35) 

The nondimensional boundary conditions are as follows. At 
the cold end of the fuel. 


n 


T(n,t) -*• T_^ 


(36) 


At the liquid surface, the conservation of mass yields 


n = 0 


r(t) = p(o) [v(o)-Xg] = 


(37) 


The conservation of energy across the surface becomes 
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(38) 


where 3 is the ratio of liquid to gas specific heats 

p 

and L is the nondimensional heat of vaporization 


(39) 


L = 
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C T* 
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( 40 ) 


The insolubility condition requires that 
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= r[0(o) - T(o) - q] + P 
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Jn=o+ 


The Clausius-Clapeyron equation after nondimensionalization is 


-L. 


0(0) = exp + T(0) 


(42) 


where 


= T = 1 

^ p* RT* 

“o 


(43) 


The dimensionless conditions far from the surface reduce to 


n oo 

0(Ti,t) T(n,t) -*• T(t) (44) 

where T(t) is a specified function* 

3. 3.3 Perturbation Analysis 

To solve the system of partial differential equations, 

Equations (23-25 and 34), with the given boundary conditions, it is 
assumed that all amplitudes of oscillation are small* This assumption 
allows the use of a perturbation analysis. For oscillatory solutions, 
the dependent variables can be expressed as a function of a perturbation 
parameter e, where e is assumed to be small 

P(t) = 1 + 

T(n,t). = T (n) + eT, (n)e^“^ 
o 1 
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o 1 
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( 45 ) 
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In these equations £ represents the normalized amplitude of the imposed 
pressure oscillations. Substituting Equation (45) into Equations 
(23-25, and 34) and the boundary conditions, the resulting equations 
are separated in terms of like powers of ee 

3.4 Zero Order (S teady State) Problem 

3.4.1 Equations 

For the zero order problem, the gas phase can be solved 
independently of the liquid phase. The problem reduces to the 
following; 


Gas Phase, 0< n < 
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with the boundary conditions 
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and at n = 0* 
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= q + T^(0) - 1 
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(48) 
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In Equations (46-49) use has been made of the conservation of mass 
requirement that 


r = constant = p v^ 
o “o “>o 


1 


(50) 
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and also that 


0^(n) - 1 (51) 

from Equations (24) and (45). 

In Equation (46), the quantity A is unknown since it is 
related to the burning rate p* which is unknown. The three 
boundary conditions permit a unique value of A to be determined. 

The liquid phase solution is required for the complete 
specification of the first order problem. To zero order in e, the 
equation and boundary conditions are: 


Liquid Phase, 


6 


1 



-oo < r\ < 0 



(52) 


n-^-co x-^T_ n = o T=T(o) 
o o o 


(53) 


The solution of Equation (52) subject to the boundary conditions 
of Equation (53) is 


T^ = T_^ + [T^(0) - T^^] exp(n/6^) (54) 

3.4.2 Solution Method 

The zero order problem is given by Equations (46-49). 
Equation (46) is highly nonlinear and must be solved numerically. 
Since the domain of interest Is half-infinite, one obviously can not 
integrate to infinity numerically. One approach to this problem is 
to integrate from T] = 0 to some finite r|, the outer flame edge, 
where T^ is within some small percentage of its true value at 
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infinity, i.e,, unity. This was the approach used in Reference 

[32], A more rigorous approach is to consider the asymptotic 

form of Equation (46) for large T]. The form of the solution for T^ 

for large n depends on the reaction order n, and further considerations 

are limited to first and second order reactions . It is shown in 

Appendix B that with n = 2, by balancing convection and reaction, 

T (n) must be of the form 
o 

A -1 

= 1 - [- exp (-E)n - C] (55) 


for large X] where C is an unknown constant. For n = 1, be 

of the form 


T 

o 


1 - C exp 


^ 1 - A+4A expT-^^ j 


(56) 


for large n where again C is an unknown constant. 

With the asymptotic form of T^ for large n. Equation (46) 
can be solved numerically without specifying an arbitrary outer flame 
edge. Equation (55) or Equation (56) is used to generate starting 
values for the numerical solution for some large, but finite, r\ with 
an assumed value of C* With these starting values, Equation (46) is 
integrated numerically to n = 0 with an assumed A. At ri = 0 the 
boundary conditions, Equations (48-49), must be satisfied and thus 
A and C are determined uniquely. However, since Equation (46) is 
nonlinear, a double iteration technique must be used and Equation 
(46) solved at each step of the iteration in A and C until Equations 
(48-49) are satisfied. With this solution technique, the outer boundary 
condition. Equation (47), is only satisfied at the true mathematical 


infinity. 
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3.5 First Order Problem 

3.5.1 Equations 

To first order in e the equation and boundary conditions 
for the liquid phase are: 


d^T dT 

,_2 ^ ” ^‘^’^1 


dT 


= r 


dn 

T, 


n = 0 


1 dTi 


T^ = T^(0) 


(57) 

(58) 


The solution of the first order system is 


T, = {T, (0) + 

1 1 o,iw o 


[T^(0) - T ^]> exp 


ri + /I + 4ici)^,* 
[ 26^ ’ 


-r 


^•[T^(0) - T_J exp(n/6^) 


(59) 


The first order gas phase equations reduce to the following. 


Gas Phase, 0 < r| < 


d^9. 


dn 


dG 


rj^ = constant 


- - i"»i “ -[^] 


iW T 


(60) 


(61) 


^1 '^^l / r(l+6-n| n E ■} 

,2 " dn ^ly^o [ T 1-T 2 ^ 

dn L' o / o T J 


P t I 

o 


dT d T .r ,-1 

^ dn^ I ^ J 


ioiT - (n-l)qw 
o ^ 0 


In Equation (62) is defined by the following expression: 


(62) 


= — T (1-T ) exp [— ] 

o n o o 1 

q o 


( 63 ) 
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Using Equation (59) the conservation of energy boundary 
condition to first order in e becomes 


n = 0 


+ T^(0)[1-B + (1 + /i+4ia)6^)] - ^ 

1 •’n= 


(0)(1-S) + L} 
(64) 


The Clausius-Clapeyron equation reduces to 




L T,(0) 

■^O 


- 1 


and the insolubility condition requires that 


(65) 


d6 ] 

^ = r Il-T (0)-q] + e (0)-T (0) + 

■ ' Jn=o+ ^ j- j. 


dT. 


dn 


n=o+ 



ri=0+ 


( 66 ) 


The outer boundary conditions from Equation (44) are 


r| 00 

^■^(n) Tj^(n) a constant (67) 

The constant in Equation (67) can be determined by solving Equation 
(61) for large q. It is easily verified that as q -»- «>, 

6i(n) Ti(n) ^ (68) 


which is the form for isentropic flow. 
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3. 5,2 Solution Method 

The first ofer problem involves the solution of two linear 
second order ordinary differential equations, Equations (61 and 62), 
involving an unknown constant, r^. Five boundary conditions, Equations 
(64-67) are imposed on the solution permitting a unique value of r^^ 
to be determined. 

The solution method to the first order problem is discussed 
in detail in Appendix C. Briefly, the method consists of exploiting 
the linearity of the system to separate the problem into the numerical 
solution of five second order ordinary differential equations. 

As in the case of the zero order problem, the equations are integrated 
to h = 0 using asymptotic analysis to generate starting values for 
the numerical routine at large ri- 

For the first order problem the constant, is related to 

the unsteady burning rate in response to the imposed pressure 
oscillations of amplitude £, However, is not measured directly 
in the experimental portion of this study. The amplitude and phase 
angle of the liquid surface oscillations, , are measured 

experimentally. However, and are uniquely related as developed 

in the following. 

From the conservation of mass at the liquid surface, 

Equation (37), and Equation (45) it follows that 

(69) 

Since is a given thermodynamic property, is directly related 

to r^. Moreover, is a complex number with both amplitude and 
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phase parameters. Both of these parameters can be directly compared 
with the measured amplitude of the liquid surface oscillations and 
the phase angle of these oscillations with respect to the pressure 
oscillations . 



CHAPTER IV 


RESULTS AND DISCUSSION 

In this chapter the results of this investigation are 
described and discussed. In the first part of the chapter the 
findings of the steady strand experimental study are presented and 
compared with the theoretical predictions. The latter part of the 
chapter deals with the experimental and theoretical results of the 
oscillatory combustion case. 

The theoretical model is described in detail in Chapter III. 
Since the experimental study was conducted with liquid hydrazine 
as the fuel, the properties used in the theoretical model are 
those applicable to hydrazine. These properties are listed in 
Table 1. The specific correlations and sources of data for these 
properties are given in Appendix D. 

4 . 1 Steady Strand Combustion 

This section deals with the steady strand results. As 
discussed in Chapter II, the steady strand experimental data 
consisted of a complete liquid temperature record displayed on an 
oscillograph, a film record of the liquid surface motion and the 
position of the thermocouple relative to the liquid surface as a 
function of time. From this data both steady strand burning rates 
and liquid phase temperature distributions were obtained. 



Table 1 


Properties Used in the Theoretical Model 




Property 

Value 


T* (K) 
—00 ' ' 

298 

6 

0 

q*(cal/gm) 

1180 

L*(cal/gm) 

410 

L*(cal/gm) 

9750 

iln a* (atm) 

12.601 

Y 

1.126 

Pj(gm/cm^) 

1.0 

(cal/gm-^K) 

0.74 

Cp(cal/gm-“K) 

0.74 

(cal/cm-sec-* K) 

coO 

4.2 X 10" 



46 


4 4 1.1 Burning Rates 

As discussed in the Introduction previous investigators 
found that the purity of the liquid fuel had a significant effect 
upon hydrazine strand burning rates. [12-15] In order to clarify 
this point, a number of tests were conducted with hydrazine of 
varying liquid purities. These results are summarized in Figure 6, 

The burning rates presented in Figure 6 were obtained using 
a 12 mm I.D# tube. While the differences between the 99,6 and 98,4% 
purity hydrazine were small, adding distilled water to the hydrazine 
as an impurity (95.6, 92,6% purity) resulted in a significant 
reduction in the burning rate in agreement with previous investigators. 
[12-15] However, Antoine [15] found that hydrazine exhibited abrupt 
increases in its burning rate for a small change in pressure. The 
present measurements did not detect any such discontinuous behavior. 

The 98-6% purity hydrazine was used in the remainder of the testing 
due to its greater availability. 

Experimental steady strand burning rate measurements were 
terminated at the upper extinction pressure. Above this pressure, 
the fuel failed to sustain combustion following ignition. The 
upper extinction pressure increased slightly as the purity of the 
fuel increased. The upper extinction pressure also increased as 
the diameter of the sample tube increased. 

Gray and Kay [13] also found lower extinction pressures 
for hydrazine. The lower extinction pressure was not investigated 
in the present study. It was observed, however, that slight bubbling 
of the liquid occurred at the lower pressures tested. Gray, et al,, 
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•-r ~T. IV. 

PRESSURE (atm) 


Figure 6 Steady Strand Burning Rates as a Function of Pressure 
for Various Liquid Purities 
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[14] indicate that bubbling of the liquid is an indication that 
the lower extinction pressure is being approached. 

It is well known that steady strand burning rates are 
dependent upon sample tube diameter due to surface tension effects 
providing more liquid surface area for the combustion process. 

This effect is shown in Figure 7 at various pressures. By plotting 
the data as a function of reciprocal tube diameter, linear fits 
were obtained, allowing a determination of the correct fundamental 
burning rate (burning rate in an infinite diameter tube) by 
extrapolation. 

Figure 8 is a plot of the fundamental burning rate as a 
function of pressure. Also shown on the plot are the theoretical 
predictions for various reaction orders and dimensionless activation 
energies. In order to obtain the theoretical curves shown on 
Figure 8 the theory was matched to the experimental results at 
atmospheric pressure assuming a first order reaction and at 6.7 
atm assuming a second order reaction. This matching was required 
since the gas phase global kinetic constants are not known for 
liquid hydrazine strand combustion. Table 2 lists the correlation 
conditions and resulting dimensional pre-exponential factors used 
in Figure 8. These pre-exponential factors were used in the 
remainder of this study, except as noted later. 

The results shown in Figure 8 indicate that the gas phase 
reaction for hydrazine strand combustion is first order at low 
pressures switching to second order at high pressures. These 
results are in agreement with the earlier work of both Gray [13, 14] 
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RECIPROCAL OF TUBE 
DIAMETER (mm'') 


Figure 7 Steady Strand Burning Rates as a Function of the Reciprocal 
of the Tube Diameter 
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PRESSURE (atm) 


Figure 8 Theoretical and Experimental Steady Strand Burning Rates as 
a Function of Pressure 



Table 2 


Correlation Conditions and Parameters Used in the Theoretical Model 


_ _ Correlation Condition „ Correlation Condition 

n - 1 n = 2 

P* = 1,0 atm; v* 0.021 cm/sec P* = 6.7 atm; v* = 0.084 cm/sec 


E 

E* 

(kcal/mole) 

B*i 

(sec 

3 

(cm /gm-sec) 

1 

2.673 

3.473 X 10^ 

7.608 X 10^ 

10 

26.73 

7.468 X 10® 

8.385 X 10^^ 

15 

40.08 

— 

3.749 X 10^^ 

30 

80.16 

— 

8.446 X 10^^ 


Ui 
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and Antoine [15], Gray found that for subatmospheric pressures the 
hydrazine strand burning rate was proportional to the square root 
of pressure, which is indicative of a first order reaction; Antoine 
found that for pure hydrazine the strand burning rate was proportional 
to pressure for pressures in the range 1-19 atm, which is indicative 
of a second order reaction. However, the results are not in agreement 
with those of Adams and Stocks [12] who found that the burning rate 
was proportional to the square root of pressure up to about 10 atm. 

As is evident in Figure 8 the gas phase activation energy 
has little effect upon the theoretical strand burning rate. This 
result is in agreement with earlier droplet decomposition flame 
studies which also found that activation energy had little effect 
on steady burning rates. [7] 

4.1.2 Liquid Temperatures 

As discussed in Chapter II, liquid temperatures of burning 
strands of hydrazine were obtained by stretching a thermocouple 
horizontally across the tube and recording its output as the liquid 
surface propagated down the tube due to combustion. In this manner 
not only liquid temperature distributions but also liquid surface 
temperatures were obtained. 

Figure 9 compares theoretical and experimental liquid 
temperatures at two pressures, 20.4 atm and 1.0 atm, as a function 
of distance from the liquid surface. The temperature readings at 
positive distances result from a liquid film formed around the 
thermocouple as it leaves the surface (due to surface tension) and 
should be disregarded. [38] 



LIQUID TEMPERATURE (®K) 



DISTANCE (mm) 

Figure 9 Theoretical and Experimental Liquid Phase Temperature 
Distributions 
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As is evident from Figure 9 , the theoretical and experimental 

results agree quite well<, The results at other pressures were similar 

to those shown in Figure 9, justifying the use of constant liquid 

phase properties • The value of the liquid thermal diffusivity used 

-3 2 

to compute the theoretical curves was 1.26 x 10 cm /sec, which is 
a reasonable value for a liquid such as hydrazine. 

The experimental values of the liquid surface temperatures 
were obtained in two ways . With the first method the surface 
temperature was taken as the temperature recorded by the thermocouple 
at the instant of time when the thermocouple was just tangent 
to the liquid surface (as obtained from the film record). With the 
second method the liquid surface temperature was taken as the 
temperature where the sharp ”knee^^ appeared in the thermocouple 
output as shown on Figure 9. The two methods yielded surface 
temperatures differing by less than 1%. 

The liquid surface temperature results are shown in Figure 10 
as a function of pressure. Two theoretical curves are also shown on 
Figure 10. The present unity Lewis number theoretical analysis 
agrees well with the data, justifying the assumption of conventional 
phase equilibrium at the surface. An infinite activation energy 
analysis is also seen to yield essentially the same results. The 
theoretical development using the infinite acitvation energy 
assumption is presented in Appendix E. 

4.1.3 Additional Theoretical Results 

Solution of the steady gas phase problem gives not only 
the nondimens ional eigenvalue. A, and the liquid surface temperature 



SURFACE TEMP. ( “K ) 



0 2 4 6 8 10 12 14 16 18 20 22 


PRESSURE (atm) 

Figure 10 Theoretical and Experimental Liquid Surface Temperatures 
as a Function of Pressure 
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but also the complete gas phase temperature and reaction rate profiles 
at a given pressure (for an assumed reaction order and activation 
energy) . In this section the theoretical gas phase temperature and 
reaction rate profiles are described. 

Figures Il-'IS show the gas phase temperature and reaction 
rate profiles for a second order reaction with nondimens ional 
activation energies of 1, 10 and 30, respectively. The nondimens ional 
reaction rate, is defined as 


= AT^ (1 - T ) exp (■ 


FO 


T„)/q 


n 


(70) 


As the activation energy increases. Figures 11-13 indicate 
that the region where gas phase reaction effects are Important 
becomes narrower and the reaction rate profiles become more sharply 
spiked. In addition the temperature reaches its limiting value 
closer to the liquid surface for a higher activation energy. 

The qualitative aspects of the first order profiles are 
similar to those of a second order reaction. First order profiles 
are shown in Figures 14-15 for nondimens ional activation energies of 
1 and 10, respectively. 

The effect of pressure on the nondimenslonal profiles for 
both a first and second order reaction was slight. Increasing 
pressure tended to cause the reaction rate profile to peak slightly 
closer to the surface with the result that the temperature more 
rapidly approached its limiting value. 
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Figure 11 Steady State Gas Phase Temperature and Reaction Rate 
Distributions for n = 2, E = 1 
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Figure 13 Steady State Gas Phase Temperature and Reaction Rate 
Distributions for n = 2, E = 30 
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Figure 14 Steady State Gas Phase Temperature and Reaction Rate 

Distributions for n = 1, E =» 1 ' S 
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Figure 15 Steady State 
Distribution 
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4.2 Oscillatory Combustion 

The results described in the previous section of this chapter 
indicate that the steady combustion portion of the theoretical model 
is in good agreement with available experimental measurements of 
the burning rates, liquid temperature distributions and surface 
temperatures of hydrazine strands. In this section experimental 
oscillatory strand burning rates are compared with the theoretical 
model. 

4.2.1 Data Reduction 

The data from the oscillatory version of the strand combustion 
apparatus consisted of an oscillograph record of the pressure in 
the test chamber and a film record of the liquid surface motion, 
both as a function of time. Typical experimental data from the 
oscillatory combustion apparatus are shown in Figure 16, 

In Figure 16 the pressure record taken from the oscillograph 
trace and the liquid surface motion data taken from the film record 
are drawn to the same time scale. These results allow the determination 
of the amplitude and also the phase angle of the surface oscillations. 

As shown in Figure 16 the phase angle of the surface 
oscillation with respect to the pressure oscillations was determined 
by measuring the distance between the peaks of the two curves. 

The amplitude of the oscillations was measured as half the distance 
between the maximum and minimum in the surface oscillations for one 
cycle. Some repeatability tests were conducted; the maximum 
difference between two sets of data taken at the same test condition 
was less than 25%. For those conditions where more than one set of 
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data was taken, the experimental points shown in the following 
figures are the average of the separate tests • 

4,2.2 Effect of Pressure Amplitude 

The theoretical model described in Chapter III assumed that 
the amplitude of the pressure oscillations was small enough so that 
only linear effects needed to be considered in the perturbation 
analysis. In order to check experimentally where this assumption 
breaks down, a series of tests were conducted with constant mean 
pressure and frequency of the imposed pressure oscillations but 
varying pressure amplitude. For these tests the mean pressure was 
held constant at 1,54 atm and the frequency was fixed at 1,32 Hz. 

In Figure 17 the amplitude of the liquid surface oscillations 
divided by the normalized pressure amplitude, e, is plotted as a 
function of pressure amplitude. The linear analysis assumes that the 
data plotted in this manner should be constant. It is evident from 
Figure 17 that the data begin to deviate significantly from this 
constant relation for pressure amplitudes above about 0.25, 

In Figure 18 the mean burning rates for the same tests as 
shown in Figure 17 are plotted as a function of normalized pressure 
amplitude. Again the linear theory assumes that the mean burning rate 
is the same for all pressure amplitudes ♦ As before, significant 
deviations from linearity appear to arise for pressure amplitudes 
above about 0*25. 

Thus, both sets of experimental results indicate that 
nonlinearities in the' combustion system are important for pressure 
amplitudes above about 0,25. To insure the validity of the linear 



Figure 17 Amplitude of the Liquid Surface Oscillations as a Function 
of Pressure Amplitude 
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Figure 18 Mean Burning Rates as a Function of Pressure Amplitude 
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theory all additional experimental data was taken with pressure 
amplitudes less than 0*15. 

4.2.3 High Pressure Results 

The main thrust of the present work was to determine the 
effect of frequency upon the oscillatory burning rate of hydrazine. 

Two sets of data were taken with variable frequency; the first 
set was taken at a mean pressure of 9.77 atm and the second set at 
a mean pressure of 1.18 atm. This section describes the results 
at 9.77 atm while the next section describes the results at 1.18 
atm. 

One of the important parameters obtained from the theoretical 
model is r. , the nondimensional mass burning rate. However, 
the dimensional amplitude and phase angle of the liquid surface 
oscillations are the quantities that are measured experimentally. 

As described in Section 3.5.2, r and x^^ are related by Equation (69). 

Figure 19 is a plot of the surface oscillation amplitude 
as a function of frequency at a mean pressure of 9.77 atm. In 
addition to the curves obtained from the theoretical model, the 
curve obtained from a totally quasi-steady analysis (neglecting 
both liquid and gas phase transient effects) is also shown on Figure 
19. The totally quasi-steady analysis is developed in Appendix E, 

The transient analysis indicates a larger amplitude of 
surface oscillation than the totally quasi-steady analysis for the 
range of frequencies shown on Figure 19. However, at lower frequencies 
the two models approach one another as they must. Note that the gas 
phase activation energy has a pronounced effect upon the predicted 
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FREQUENCY (Hz) 


Figure 19 Theoretical and Experimental Magnitude of the Liquid Surface 
Oscillations as a Function of Frequency at a Mean Pressure 
of 9.77 atm 
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transient results* The results for a dimensionless activation 
energy of 30 do not differ appreciably from the results for E = 15 
in this frequency range and these two curves are shown as one, labelled 
E = 15, As shown in Figure 19, the theoretical curve with a 
dimensionless activation energy of 15 is in reasonably good agreement 
with the experimental results* 

Figure 20 is a plot of the phase angle of the surface 
oscillation as a function of frequency at a mean pressure of 9*77 
atm. At low frequencies, both the transient and the quasi-steady 
analyses indicate a 90° phase angle* At this condition the burning 
rate and pressure oscillations are in phase. As seen from Equation 
(69), for the burning rate in phase with the pressure, the surface 
oscillation must be 90° out of phase* Of course, the quasi-steady 
model does not predict any deviation from the 90° out of phase 
condition as frequency increases* 

As in the case of the amplitude of the surface oscillation, 
the activation energy has a significant effect upon the predicted 
phase angle. Again a nondimens ional activation energy of 15 agrees 
with the data reasonably well as shown in Figure 20. Note that the 
phase angles at these conditions are greater than 90° indicating that 
the burning rate oscillation was leading the pressure oscillation. 

Since the steady state results, as shown in Figure 8, 
indicate that the gas phase reaction at 9.77 atm is second order, 
no theoretical predictions were made assuming a first order 
reaction at this mean pressure* A second set of data was taken at 1.18 
atm where, from Figure 8, the gas phase reaction appears to be first 


order. These results are described in the next section. 
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Figure 20 Theoretical and Experimental Phase Angle of the Liquid Surface 
Oscillations as a Function of Frequency at a Mean Pressure 
of 9.77 atm 
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4.2.4 Low Pressure Results 

In order to investigate the effect of mean pressure on the 
oscillatory burning rate of hydrazine, a second series of data were 
taken at a mean pressure of 1.18 atm. At this mean pressure the 
steady results, as shown in Figure 8, indicate that the gas phase 
reaction is probably first order. 

In Figure 21 the amplitude of the liquid surface oscillation 
is plotted as a function of frequency. The theoretical curves 
obtained from the transient model assume a first order reaction. 

At these conditions, the theoretical curve for E = 1, predicts a 
smaller surface oscillation amplitude than the totally quasi-steady 
analysis. However, both theoretical curves approach the quasi-steady 
results for low frequency. 

The results shown in Figure 21 suggest that a higher activation 
energy than E = 10 is required to agree with the experimental data. 
However, the opposite conclusion is reached if the phase angle is 
considered as shown in Figure 22. The phase angle results suggest 
that a low activation energy, between 1 and 10, would best agree 
with the experimental data. Thus, the first order reaction 
theoretical results are not in good agreement with the experimental 
data at a mean pressure of 1.18 atm. 

If the reaction is assumed to be second order at 1.18 atm, 
the theoretical steady strand burning rate is significantly less 
than the experimental result as seen from Figure 8. To adjust the 
steady burning rate, the results shown in Figure 8 were 
recorrelated by matching theory to the experimental results at 1.0 



Figure 21 Theoretical (n = 1) and Experimental Magnitude of the 

Liquid Surface Oscillations as a Function of Frequency at 
a Mean Pressure of 1.18 atm 
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Figure 22 Theoretical (n = 1) and Experimental Phase Angle of the 

Liquid Surface Oscillations as a Function of Frequency at 
a Mean Pressure of 1.18 atm 
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atm for a second order gas phase reaction rather than at 6.7 atra. 

The parameters obtained from this match are listed in Table 3. 

In Figure 23, the amplitude of the liquid surface oscillation 
is plotted as a function of frequency. A second order gas phase 
reaction with an activation energy of between 10 and 15 agrees well 
with the experimental data. The same conclusion is reached when 
considering phase angles as shown in Figure 24, Thus, as with the 
high pressure results, the theoretical model assuming a second order 
reaction with a nondimens ional activation energy of about 15 is in 
good agreement with the experimental data. The corresponding 
dimensional activation energy of about 40 kcal is also in good 
agreement with the apparent activation energy of 36 kcal reported by 
McHale, et al., [49] from a kinetic study on hydrazine decomposition. 

4.2.5 Additional Implications of the Theoretical Model 

The preceding discussion has indicated that the theoretical 
model is in good agreement with all available experimental results 
on the oscillatory burning of hydrazine. After gaining this 
confidence in the model, the model was used to predict the effects 
of various parameters of interest to engine designers on the 
oscillatory behavior of hydrazine. This section describes these 
reults in detail* 

4. 2.5.1 Effect of Mean Pressure 

One of the prime quantities of interest in the unsteady 
combustion case is the dimensionless oscillatory mass burning rate 
perturbation, due to an oscillatory pressure perturbation of 
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Table 3 


Additional Correlation Condition and Parameters 
Used in the Theoretical Model 


£ 

E* 

(kcal/mole) 

n = 2 

Correlation 

Condition 

= 1*0 atm 
o 

= 0.021 cm/sec 


B* 

(cm3/ gm-sec) 


1 

2.673 


1.806 X 10® 


10 

26.73 


2.353 X 10^^ 


15 

40.08 


1.052 X 10^^ 



24 


30 


80.16 


2.375 X 10 


AMPLITUDE OF SURFACE OSCILLATIONS 
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Figure 23 Theoretical (n = 2) and Experimental Magnitude of the 

Liquid Surface Oscillations as a Function of Frequency at 
a Mean Pressure of 1.18 atm 
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Figure 24 Theoretical (n = 2) and Experimental Phase Angle of the 

Liquid Surface Oscillations as a Function of Frequency at 
a Mean Pressure of 1.18 atm 
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amplitude e. The next nine figures indicate the effect of mean 
pressure on the behavior of for various activation energies and 
a second order reaction* 

Figures 25-27 show the amplitude and phase angle (with 

respect to the pressure oscillation) of r^ for a dimensionless 

activation energy of 1 and mean pressures of 1, 10 and 100 atm, 

respectively, as a function of the dimensionless frequency of the 

pressure oscillation. At low frequencies the burning rate oscillation 

is in phase with the pressure oscillation, and the entire process, 

both liquid and gas phase, is quasi-steady. As the frequency 

approaches the characteristics frequency of the liquid phase thermal 

wave both the phase (with respect to the pressure oscillation) and 

the amplitude of r^ begin to deviate from their quasi-steady value. 

The nondimensional characteristic frequency of the liquid phase 

-4 -2 

thermal wave increased from about 10 at 1 atm to about 10 at 100 
atm. 

As the frequency increases beyond the characteristic 
liquid thermal wave frequency, gas phase transient effects become 
important. On Figures 25-27 two theoretical curves are shown. The 
solid curve neglects all transient gas phase effects but includes 
the effects of a transient liquid phase. The dashed curve includes 
the effects of both transient gas and liquid phases. As indicated 
in Figures 25-27 gas phase transient effects only become important 
for nondimensional frequencies on the order of 0.10. The quasi-steady 
gas phase analysis was obtained by taking the formal limit O) 0 in the 
gas phase as described in Appendix C. 
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Figure 25 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 1 atm, n = 2, E = 1 
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Figure 26 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 10 atm, n = 2, E = 1 


PHASE ANGLE OF r, 




MAGNITUDE OF 



Figure 27 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 100 atm, n = 2, E = 1 
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For E = 1, Figures 25-27 indicate that as the pressure 
increases the peak in the response curve due to liquid thermal 
effects becomes more pronounced. At low pressures the phase lag is 
negative for all frequencies of oscillation. As the pressure 
increases, however, the phase lag is at first positive and then 
goes negative as o) increases. This effect was found to be due to 
the shape of the hydrazine vapor pressure curve. Theoretically, 
the liquid surface temperature oscillation amplitude is relatively 
small o Thus, at low pressures a finite change in the pressure 
causes a relatively small change in the fuel mass fraction at the 
surface. At high pressures, the same change in pressure causes a 
significant change in the surface fuel mass fraction. 

Similar trends were found for the case of E = 10 as shown 
in Figures 28-30. However, in this case the peak in the response 
curve becomes very pronounced as pressure increases. 

For the case E - 30 somewhat different behavior was found. 

As illustrated in Figures 31-33 the peak in the response curve 
becomes less pronounced as pressure increases. However, as before, 
the phase angle tends to positive values as the pressure increases 
for a given frequency. 

4. 2. 5. 2 Effect of Conditioning Temperature 

Another parameter which could be of importance with regard 
to oscillatory combustion response is the fuel conditioning temperature 
T’^^. The effect of increasing from 298 to 350 K was investigated 
at a mean pressure of 1 atm for a second order reaction with a 
nondimensional activation energy of 10. The results of these 
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Figure 28 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 1 atm, n = 2, E - 10 
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Figure 29 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 10 atm, n = 2, E = 10 
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Figure 30 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 100 atm, n = 2, E = 10 
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Figure 31 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 1 atm, n = 2, E = 30 
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Figure 32 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 10 atm, n = 2, E = 30 
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Figure 33 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 100 atm, n = 2, E = 30 
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calculations are shown in Figure 34. Figure 34 was obtained by fixing 
all dimensional quantities but increasing the fuel conditioning 
temperature by 52 K, Of course, in order to satisfy overall 
conservation of energy the flame temperature had to be increased a 
corresponding amount. Increasing the flame temperature decreased 
the values of some of the nondimens ional parameters of the problem, 
including the nondimens ional activation energy. 

By comparing Figure 34 with Figure 28 it is seen that the 
effect of increasing T ^ is to decrease somewhat the first peak in the 
r^ response curve. However, the net overall effect is rather small. 

4. 2. 5. 3 Acoustic Admittance 

The major parameter of interest to the designers in 
determining the stability characteristics of a combustion chamber is 
the real part of the acoustic admittance of the burning fuel. This 
quantity can be obtained from the present analysis as described in 
Appendix G. 

The acoustic admittance varies with the distance from the 
liquid surface, but approaches a constant value for very large 
distances. As developed in Appendix G, the value shown here is the 
value at the true mathematical infinity* 

Physically, the real part of the acoustic admittance 
represents that part of the gas phase velocity fluctuation which is 
in phase with the pressure fluctuation. Figures 35-37 show the 
real part of the acoustic admittance plotted as a function of frequency 
for various pressures for a second order reaction and nondimens ional 
activation energies of 1, 10 and 30, respectively. 
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Figure 34 Burning Rate Oscillation Amplitude and Phase Angle as a 
Function of Frequency for P* = 1 atm, n = 2, E = 9.628, 
T* = aso^K ° 
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Figure 35 Real Part 
Frequency 
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As shown in Figure 35, the real part of the acoustic admittance 
is zero for low frequencies where the combustion process is quasi- 
steady* In this limit the mass burning rate for a second order 
reaction is linearly proportional to pressure. However, from the 
perfect gas relation, the gas phase density is also linearly 
proportional to pressure; therefore, the gas phase velocity is 
independent of pressure. Thus, the gas phase velocity fluctuation 
due to a pressure fluctuation is zero at low frequencies* 

As the frequency increases, liquid phase thermal wave 
effects become important causing the real part of the acoustic 
admittance to deviate from its quasi-steady value. The curves shown 
in Figure 35 reach a maximum near the characteristic frequency of 
the liquid phase thermal wave then become negative as gas phase 
transient effects become important. The curves reach a minimum 
then become positive again for very high frequencies. 

Reference [16] indicates that for instability to occur the 
real part of the acoustic admittance should be on the order of one. 

As shown in Figure 35 the real part of the acoustic admittance for 
E = 1 only becomes of the order one at very large frequencies for 
the pressures shown. Liquid phase thermal wave effects are not of 
sufficient magnitude to cause instability for E == 1. 

Figure 36 shows a similar plot for E = 10. This figure 
indicates that as pressure increases a frequency range is reached 
where the real part of the acoustic admittance is one or greater. 

Also, with increasing pressure, the peaks in the curves due to 
liquid thermal wave effects move to higher frequencies. 
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Somewhat different behavior is evident from Figure 37. 
For E = 30, the real part of the acoustic admittance is of the 
order one or greater only for low pressures except for very high 
frequencies . 



CHAPTER V 


SUMMARY AND CONCLUSIONS 


5 • 1 Summary 

The present study was concerned with the strand combustion 
of hydrazine under both steady and oscillatory conditions. The 
specific objectives of the study were as follows: 

1. Investigate the steady strand combustion characteristics 
of hydrazine experimentally to determine steady strand 
burning rates as functions of pressure, liquid 
concentration and tube diameter. Also investigate the 
liquid temperature distributions and liquid surface 
temperatures as a function of pressure during steady 
strand combustion. 

2. Study the oscillatory combustion characteristics of 
hydrazine by experimentally determining the oscillatory 
burning rate as a function of mean pressure, and the 
frequency and amplitude of the imposed pressure 
oscillations . 

3. Develop a theoretical model of the strand combustion 
system capable of predicting the observed experimental 
results under both steady and oscillatory conditions. 

For the experimental program a steady strand combustor 
was used in the initial phases of the investigation. With this 
apparatus steady strand burning rates were measured over 
the pressure range of 0.32-42 atm and liquid temperatures were 
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obtained in the range of 1.0-20.5 atm. Liquid concentrations of 
from 92,6”99.4% were studied in tube sizes of 4, 8 and 12 mm* 

A modified version of the steady strand apparatus was 
constructed to provide a means of studying the oscillatory combustion 
of hydrazine. With this apparatus oscillatory burning rates were 
obtained in the mean pressure range of 1-10 atm with frequencies 
varying from 0.4-5. 2 Hz and pressure amplitudes up to 35 % of the 
mean pressure. 

A theoretical model of the system was developed. The major 
assumptions in the theoretical analysis were that the system was one 
dimensional, a one-step, irreversible reaction occurred in the gas 
phase and conventional low pressure phase equilibrium was adequate 
to describe the state of the liquid surface. While the liquid 
phase was assumed to have constant physical properties, a variable 
property gas phase solution was obtained through the use of a 
modified Howarth transformation. The resulting equations describing 
the combustion system were solved through the use of a perturbation 
analysis with the amplitude of the imposed pressure oscillations 
taken as the small perturbation parameter. 

The zero order model described the steady strand combustion 
system. The predictions of this model were compared to the 
experimental steady strand results. The excellent agreement between 
theory and experiment justified the extension of this theory to 
include unsteady effects. 
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5*2 Conclusions 

5.2,1 Steady Strand Combustion 

The major conclusions of the steady strand portion of the 
study are as follows: 

1, Addition of water to the liquid hydrazine reduces the 
burning rate. However, liquid purities of 98.6 and 
99.4% yielded essentially the same burning rates. 

2. The effect of tube diameter scales as the reciprocal 
of the tube diameter. By plotting steady strand 
burning rates as a function of the reciprocal of the 
tube diameter, linear plots were obtained allowing 
the determination of the fundamental strand burning 
rate. 

3, The experimental results indicated that the strand 
burning rate has a pressure dependence of 1/2 for 

sub atmospheric pressures and a pressure dependence of 1 
for pressures above atmospheric. This finding is in 
agreement with the results of previous investigations, 
[13-15] The theoretical model was matched to these 
results by assuming a first order decomposition reaction 
for sub atmospheric pressures and a second order 
decomposition reaction for pressures greater than 
atmospheric. 

4. The gas phase activation energy has very little 
effect upon the theoretically predicted steady strand 
burning rate for both a first and second order reaction. 
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5* The theoretical model, assuming constant liquid phase 
properties, agrees well with the experimental liquid 
temperature distributions. 

6. The theoretical model is also in excellent agreement 
with the experimentally determined liquid surface 
temperatures justifying the low pressure phase 
equilibrium assumption. 

5.2.2 Oscillatory Strand Combustion 

The major conclusions of the oscillatory strand portion of 
this investigation are as follows: 

1. As opposed to steady combustion, for oscillatory 
combustion the gas phase activation energy has a 
significant influence upon the predicted theoretical 
oscillatory burning rate. The theoretical results 
concerning both the amplitude and the phase angle 
of the liquid surface oscillations are in good 
agreement with the experimental results when a second 
order reaction with an activation energy of about 

40 kcal (E = 15) is assumed. 

2. Based upon the theoretical acoustic admittance, the model 
predicts that for activation energies in the range of 
25-40 kcal (E = 10-15), the hydrazine combustion system 
exhibits a range of frequencies where instability 

could occur. As pressure increases, the instability 
range moves to higher frequencies and the system becomes 


less stable. 
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A, 1 Estimation of Radiation Effects During Strand Burning 

An estimation of the effect of radiation can be made by 
comparing the energy reaching the liquid surface due to radiation to 
the total amount of energy reaching the liquid surface. In the steady 
state the total amount of energy reaching the liquid surface is equal 
to the product of the liquid density, the fuel burning rate and the 
heat required to raise the temperature of the liquid to the surface 
temperature and vaporize the fuel at the surface temperature. 

To get a rough estimation of the amount of heat reaching the 
liquid surface due to radiation, the combustion system was taken to be 
an infinite cylinder of the decomposition flame products at the 
adiabatic flame temperature of the fuel. The walls of the cylinder 
were assumed to be transparent. The shape factor and mean beam 
length were taken from Hot tel and Sarofim [40] . 

The fuel decomposition products were taken to be those 
suggested by Audrieth and Ogg [39]; namely 

^ NH^ + I ^ H 2 (A.l) 

The only emitter is the system in with a mole fraction of 0.5. 
Using the corresponding partial pressure and the mean beam length, the 
emissivity was taken from Hottel and Sarofim [40]. The heat due to 
radiation reaching the liquid surface was then the product of the 
shape factor, the gas emissivity, the Stefan-Boltzmann constant and 
the adiabatic flame temperature raised to the fourth power. 

By comparing the radiation heat transfer to the total heat 
transfer in the above manner, it was found that the radiation effect 


was always less than 1% of the total. This result is 
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not surprising since hydrazine has a low adiabatic flame 
temperature . 


A, 2 Phase Equilibrium Assumption 

The equilibrium assumption at the liquid surface can be 
justified if the rate at which molecules from the gas phase strike 
the liquid surface is much greater than the fuel burning rate* The 
number of molecules striking the liquid surface per unit time was 
estimated using kinetic theory. 

The rate at which molecules from the gas phase strike the 
liquid surface was estimated using the following equation derived 
from kinetic theory [41] 


0 = 


' 2ttRT 


(A. 2) 


where 4 is the rate at which molecules strike the surface per unit 
area, P is the pressure, M is the molecular weight and T is the 
surface temperature* 

Using the experimentally determined surface temperature as 
a function of pressure in Equation (A.2), $ was compared to the 
experimental burning rate at various pressures* It was found that 0 
was always between two and three orders of magnitude greater than 
the burning rate at any given pressure* Thus the assumption of 
phase equilibrium at the surface is justified* 


A* 3 Constant Total Pressure Assumption 

The constant total pressure assumption is really two 
assumptions, namely that the Mach number of the flow is much less 
than unity and that the wavelength of any periodic pressure disturbance 
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is long compared to the distance from the outer edge of the reaction 
zone to the liquid surface. The Mach number was estimated to be in 
the range of 0,001-0.002 and is much less than unity. This section 
of Appendix A, therefore, deals with estimating the frequency where 
the long wavelength assumption breaks down. 

To estimate this frequency assume that the temperature 
profile is linear between the surface temperature and the flame 
temperature. Then the total heat transfer to the liquid surface can 
be expressed as 




p* v| L* 


(A.3) 


All quantities in Equation (A.3) are known either as thermodynamic 

properties or as experimentally measured parameters except the distance 

from the outer edge of the reaction zone to the surface, (x^ - x*) . 

r s 

The parameter, (x| - computed from Equation (A.3) 

as a function of pressure. It was found that (x| - x^) is greater 
than 0.10 of the wavelength of any periodic pressure disturbance 
for frequencies above about 50,000 Hz, 



APPENDIX B 


FORM OF THE ZERO ORDER ENERGY EQUATION FOR LARGE n 
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The only equation which must be integrated for the zero 
order problem is Equation (46) which is rewritten here. 


d T dT , 

o o , . 1-n „ 1+o-n,, „ sn r-E, ^ 

dn ■ o 


To find the form of Equation (B.l) for large ri let 


T = 1 - ?T 
o o 


( 1 ) 


(B.2) 


where ^ is a small parameter. Combining Equation (B.l) and Equation 
(B.2) yields 

d^T dT . 

C - S — §r Aq^ ” exp(-E)T = 0 (B.3) 

dn 

The solution of Equation (B*3) depends on the order of the 
reaction, n. For n - 1, the solution of Equation (B*3) is 
straightforward . 

nj (B.4) 

The growing exponential part of the solution to Equation 
(B*3) is rejected since must approach zero for large ri from 

Equation (47). in Equation (B.4) is an unknown constant. 

Combining Equation (B.2) and Equation (B.4) results in 

- 1 - 0 nj (B.5> 

For n = 1, where 


C = € 


(B.6) 


and is an unknown constant. 
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For n = 2, Equation (B.3) can be solved by stretching the X] 
coordinate such that 


n = ^ y 


(B.7) 


where a is a constant determined by proper matching- 

Substituting Equation (B.7) into Equation (B,3) yields 

2 ( 1 ) ( 1 ) 

^l-2a o^ _^l-a _o _ ^2^ exp(-E)(T /q = 0 (B.8) 

dy ° 

A proper match of the terms in Equation (B.8) is only obtained by 

balancing convection and reaction; the diffusion term is of higher 

order. Then a = -1 and the resulting equation is 

dT 2 

Aexp(-E)(T/«) (B.9) 

solving Equation (B.9) 

T = (7exp(-E)y - C^)"^ (B.IO) 

o q 1 

where is an unknown constant. 

From Equation (B.2) and Equation (B.7) 

= 1 - ?/(^ exp(-E)n - C^) (B.ll) 

or 

= 1 - (| exp(-E)n - C)"^ (B.12) 

for n = 2, where 


C - Cj/5 


(B.13) 


and is an unknown constant. 
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The first order problem Is given by Equations (60-66) and 
Equation (68) ♦ The problem can be simplified somewhat by noting that 
the solution of Equation (61) can be represented by 

exp[ ^ + e^p(n) (c.i) 

where is the unknoxm constant in the homogeneous solution and 
represents the particular solution. In arriving at Equation (C.I) the 
growing exponential part of the homogeneous solution was rejected 
since, physically, neither the temperature nor the fuel mass fraction 
fluctuations can grow exponentially with increasing T), In general 
0^p must be determined numerically. However, for large n, using 
Equation (47) in Equation (61), 

©ipCn) ^ ^ (C.2) 

Equation (C*2) also follows from Equation (68). 

Using Equation (C.I) in Equation (62) the first order 
problem reduces to the following equation 

d^T dT- 

— ^ + pTi = r^b + K^c + d (C.3) 

dn 

where 


p = qw 


h+6-n 

" + ® 1 

T , 

Li 0 i 

1-T 2 

' ^ o T 


- itb 


b = 


2 

d T 


nqw^ 

c = - exp 

o 


dT 



dr) 

fl - v^l+4~ico 


nqw 


^ - -p 

dr) ' 


T ®1P 


(C.4) 

(C.5) 

(C.6) 

(C.7) 
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The boundary conditions can also be reduced to the following: 




(C.8) 


Ki + Qj^pCO) = jT^(O) + k 

h* " “'i ^ l^ln-o+ 


(G.9) 


where 


f = ^ [T (0) - T ] [i/l+4ia)6. - 1] + T (0)(l-e) + L (C.ll) 

26 0 1 o 

g = 1-B + (1 +A+4iw6^) (C.12) 


h = - 


(C.13) 


aqL- -L- 

j = 2 exp 1 

T^^(O) o''^'^ 

k = T (0) - 1 
o 


(P.14) 


(C.15) 


£ = (1 - /i+4iw) /2 


(C.16) 


m = 1 - T (0) - q 
0 


(C.17) 


The outer boundary condition remains as specified by Equation (67) 
or Equation (68) . 

The equations representing the quasi-steady gas phase 
problem can be obtained by setting O) *= 0 in Equations (61 and 62). 
These equations can be summarized as follows: 


®i"'i 


(C.18) 



since 


01P = 0 


and 


d dT^ 


+ p T, = r,b + K,c + d^ 


dn 


2 dn si 1 Is s 


(C.19) 


where 


p = qw 
s o 


d = - 
s 


d^T 


dn 



(C.20) 

^■^o T 


nqw^ 0 


l-T 

(C.21) 

0 


- qw^(n-l) 

(C.22) 


The boundary conditions remain the same except 

/de, ‘ 


I dn 


IP 


= 0 


'n=o+ 


(C.23) 


Since the equations for the complete unsteady and quasi-steady 
problems are similar except for the values of the coefficients, the 
general solution method as discussed in the following is applicable 
to either case. 

The first order problem reduces to the solution of one 
linear, second order, ordinary differential equation. Equation (C.3), 
with two unknown constants r^^ and The four boundary conditions 

permit unique values of the unknown constants to be found as well 
as the desired solution out of the family of solutions of the 
second order, ordinary differential equation. 

Exploiting the linearity of Equation (C.3), let 
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Then Equation (C.3) can be separated into three equations 


d dT^ 


dn^ 

dn 

d^T* 

dT* 


" dn 

2- 


dTi 

dTi 

dn^ 

dn 


+ pT^ = d 


+ pT* = b 


+ pT^ = c 


The boundary conditions, Equations (C.8-C.10 and 67) can also be 
separated as follows; 





gT^(O) + h 


dT* 


n=o+ 


= gx*(0) + f 


dTj 


n=o+ 


= gX^(O) 


e^p(o) = jT^(o) + k 


0 = jX*(0) 


1 = jX^(O) 



6ip(0) - X^(0) + 


y\ 



0 = m - T*(0) + 



(C.25) 

(C.26) 

(C.27) 

(C.28) 

(C.29) 

(C.30) 

(C.31) 

(C.32) 

(C.33) 

(C.34) 


n=o+ 


(C.35) 
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and as r| -»■ “ 


/dT 

t . 1 . 1^(0) + 


Tl-C 


'tt=0+ 


T* ^ C* 


-f C 


(C.36) 


(C.37) 

(C.38) 

(C.39) 


where C, C* and C are constants. 

' The three differential equations. Equations (C.25-C.27), 
each have a particular and a homogeneous solution. However, since 
the homogeneous equations are identical, the homogeneous solutions 
are Identical. Breaking the equations into homogeneous and particular 
parts, let 


Tip * 

(C.40) 

T*p + 

(C.41) 


CC.42) 


where K, K* and K are unspecified constants. The differential 
equations, Equations (C.25-C-27) become 


d^T dT 


dn 


/V A 

-^ — -P-1P = ^ 


(C.43) 


(C.44) 
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d dT* 

IP IP 

±£ ±L 4* = b 

dn' 


(C.A5) 


d^T 


IP 


dT 


IP 


dn 


dn 


+ pT 


IP 


(C.46) 


The constants K, K* and K are chosen such that one of the boundary 
conditions at ri=0, Equations (C,8-C.10) is satisfied. The remaining 
two boundary conditions at = 0 are used to determine r^^ and 

At this point one must choose which boundary condition to 
use to determine K, K* and K. The solution of the system of equations, 
Equations (C. 43-46) should be unique and thus independent of the choice 
of boundary condition ♦ It was verified that the choice has no 
effect on the final computed results by using all three boundary 
conditions, Equations (C.8-C-10), in turn and computing, numerically, 
the solution. Equation (C.9) is chosen as the boundary condition 
to be used here for illustration purposes. 

Separating Equations (C.31-C.33) using Equations (C,40“C,42) 


yields 




. 8j^p(0) - k - jTj^p(O) 

K* = - T^(0)/T^^(0) 

1 - jT (0) 


(C.47) 

(C.48) 

(C.49) 


Summarizing thus far, the solution method consists of 

✓N 

solving Equations (C,43-C.46) numerically for T*p(r)) 
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and T^p(ri). Equations (C.47-C-49) are then used to determine K, K* 
and K. Next, Equations (C.40-C.42) are used to find T^(ri), T*(ri), 
and T^(ri)* Finally, Equations (C.28-C.30) and Equations (C.34-C.36) 
are recombined to produce two equations involving two unknown 
constants, and These equations are: 


A 



(C,51) 


Equations (C.50-C*51) permit unique values of r^ and to be 
determined. 

The first order problem is reduced to the solution of five 
second order, ordinary differential equations, Equations (C.43-C.46) 
and the solution of Equation (61) for with the zero order 

problem, these equations were solved by integrating to the liquid 
surface from some large, but finite, n. Asymptotic analysis was 
used to determine the analytical form of the solutions for large n* 
These analytical equations were then used to generate starting values 
for the numerical solution. The asymptotic analysis is discussed 
in the following. 

The form of 0j^p(Ti) for large n has already been discussed. 
For the complete unsteady gas phase problem, the form is given by 
Equation (C.2). For the quasi--steady gas phase problem the form of 
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for all n is given by 


©IP = 0 (C.52) 

For Tj^p, T*p and f^p the form of the solutions depends 

on the reaction order, n. For n = 1 and 6 = 0, it can be shown that 
for large n 


-A exp(-E) - iW = -p^ 

(C.53) 

-(-^)iti) - Aexp(-E)0^p 

(C.54) 

b 0 

(C.55) 

-Aexp(-E) exp[-^^ /l+4iw 

(C.56) 


Substituting Equations (C.53-C.56) into Equations (C.43-C.46) it is 
found that for large n 


^IH “ 


1 - v^l+^ip' 


n] 


= exp[- 


1 - V14?F 


n] + 


+ Aexp(-E)6 


IP 


Aexp(-E) + iw 


(C.57) 

(C.58) 


T* = T 
IP IH 


(C.59) 


T = 
IP 


1 - v^+4p^ 
exp[ r n] + 


piexp(-E) + 

1 - /l+4i(iJ 
2 

2 

1 - v^l+4iU 
2 

Aexp(-E) + ±0) 


fl - ✓'1+4 iW 
expt n n] 


(C.60) 


For n = 2 and 6 = 0, it can shown that for large r| 


p ^ - ico 


(C.61) 
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In deriving Equations (C.57-C.60) and Equations (C.65-C.68) 
growing exponentials were disregarded in order to satisfy Equations 
(C.37-C.39). The constants in the homogeneous solutions were 
normalized to unity since they were already considered in the constants 
K* and K. 

Summarizing, Equation (61) is first solved to find using 

Equation (C.2) to generate starting values for large n* If the gas 
phase is assumed to be quasi-steady then is set equal to zero for 
all r|. Next, Equations (C.43-C.46) are solved using the appropriate 
starting values. Equations (C.57-C.60) for n = 1 or Equations 
(C.65-C.68) for n = 2. Equations (C.47-C.49) are then used to find 

A. /\ 

K, K* and K. Next T^, T*, and are found using Equations (C,40-C.42). 
It should be mentioned at this point that the first and second 
derivatives of T^, T* and can be found from equations similar to 
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Equations (C.40-C.42). and are then determined from Equations 

(C.50-C.51). Finally, Equation (C.24) is used to find Tj^. 



APPENDIX D 


PHYSICAL PROPERTIES 
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The references for the physical properties required in 
the calculations are shown on Table 4* The correlations used to 
compute the properties are discussed in the following. 

The specific heat of the fuel in both the gas and liquid 
phases were assumed constant. The gas phase specific heat was 
evaluated at a temperature of 900*^K. 

The gas phase thermal conductivity was evaluated at 
the adiabatic flame temperature. For this calculation the gas 
phase mixture was assumed to be made up of 50% ^ 2 ^ 4 ^ 12.5% 

and 12.5% on a molar basis. 

The liquid phase thermal conductivity was assumed constant. 

It was evaluated from the liquid thermal diffusivity which was 

found to give the best fit of the zero order theory to the experimental 

results on liquid phase temperature distributions. This value of 

-3 2 

liquid thermal diffusivity was found to be 1.26 x 10 cm /sec 
as discussed in Section 4.1.1. 

The adiabatic flame temperature, was calculated allowing 

for all relevant dissociation reactions. The thermo chemical 
properties required for this calculation were taken from the JANAF 
Tables. [46] 

The heat of vaporization of the fuel was evaluated from the 
overall conservation of energy relation. By specifying the temperature 
at the cold end of the liquid, the final gas phase temperature, the 
heat capacities of the gas and liquid, and the heat of combustion, 
the overall conservation of energy in the steady state could only be 
satisfied by proper selection of the heat of vaporization. The value 
used is shown in Table 1. 
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Table 4 


References for Physical Properties 





References 


Property 

«2«4 

Combustion Products 

h° 

42 


43 


39 


— 

X* 

TOO 

44^ 


44 

C 

P 

45 


— 


39 


— 


^Computed, Method of Reference [44] 
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The constants in the Clausius-Clapeyron equation were 
determined from a least squares fit to the existing data on vapor 
pressures as a function of liquid temperature for hydrazine, 

[45, 47-48] The values found are listed in Table 1. 



APPENDIX E 


SURFACE MASS FRACTION USING THE INFINITE 
ACTIVATION ENERGY ASSUMPTION 
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If an infinite gas phase activation energy is assumed rather 
than a unity Lewis number, the zero order equation relating surface 
fuel mass fraction to surface temperature, Equation (16), is 
different. In this Appendix, the theoretical development of an 
equation analogous to Equation (16) is presented for an infinite 
activation energy flame. 

The theoretical development follows closely that given in 
Reference [38], Faeth [38] used a molar basis for his analysis 
and assumed that both the specific heat and thermal conductivity of 
the gas phase were linear functions of temperature. In the 
analysis presented here, a mass basis is used and the specific heat 
is assumed constant but thermal conductivity is a function of 
temperature . 

Since the reaction is confined to a flame surface, the 
steady state energy and species conservation equations can be written 
as 


p*v*[C (T^ - T*) + L*I - A* ^ (E.l) 

P S 

P*v*[Yp - 1] = p*D ^ (E.2) 

With the assumption of infinite activation energy, the outer 
boundary condition can rigorously be applied at the outer edge of the 
flame zone, x|, since the flame is infinitely thin. 

Thus , the boundary conditions are 


x*=0 Y=Y T=T* 

F ^F,S s 


Y 


F 


T = T* 
00 


X* = X* , 


0 


(E.3) 
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The quantity X*/p*D is only a weak function of composition 
and is assumed constant. [38] 

Dividing Equation (E.l) by Equation (E.2) to eliminate spatial 
derivatives and integrating, using the boundary conditions, Equation 
(E.3) yields 

X*/(p*DC ) 

Y = 1 _ r L* ' P 

F,S (T*-T*) + L*^ (E.4) 

p s 

Equation (E.4) then replaces Equation (16) if an infinite activation 
energy is assumed. 




APPENDIX F 


DEVELOPMENT OF THE TOTALLY QUASI-STEADY ANALYSIS 
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The totally quasi-steady analysis neglects both transient 
liquid and gas phase effects. The analysis was done by using the 
experimental data on the fundamental strand burning rate of hydrazine 
as a function of pressure shown in Figure 8 and fitting it to the 
following empirical equation for pressures greater than atmospheric. 

dx* 

V* . u» + uj PJ - <F.l) 

where and u* are constants, 
o 1 

Assume the pressure^ P*, varies sinusoidally with time such 

that 

p* = p;fe -h p* <F.2) 

o oo ol 

where and are also constants. Substituting Equation (F.2) into 
(F.l) and integrating yields 

X* = (u* + u* P* )t* + sin w*t* (F.3) 

Equation (F.3) is the totally quasi-steady theoretical prediction of 
liquid surface movement. The first term on the right hand side of 
Equation (F.3) represents the steady regression of the liquid surface 
for constant pressure. The second term represents the oscillatory 
motion of the liquid surface under oscillatory pressure conditions. 

Note that the amplitude of the oscillating component of the liquid 
surface motion is an inverse function of frequency, and that the 
phase angle of the surface oscillation with respect to the pressure 
oscillation is constant and equal to 90°. 



APPENDIX G 


DEVELOPMENT OF THE ACOUSTIC ADMITTANCE EXPRESSION 
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In order to determine the acoustic admittance of the combustion 
system, the gas phase velocity perturbation due to a pressure 
perturbation of amplitude e must be known* The solution of the 
equations developed from the present theoretical model do not explicitly 
yield the gas phase velocity perturbation. However, this parameter 
can be found from other output parameters of the theory as developed 
in the following. 

Rewriting Equation (21) here for convenience 

X 


r 



I 


X^(t) 


p dx 


(G.l) 


but 



(G.2) 


and therefore 


rT = Pv + Tj dx- (G.3) 

Perturbing Equation (G.3) using Equation (45) and rearranging 
yields to zero order in E 


V = T 

o 


(G.4) 


and to first order in e 
where to first order 



iwx ,T 
si Q 

T^(0) 


(G.5) 


dx = T dn 
o 


(G.6) 


and 
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n 

$ = dr] 


The acoustic admittance is defined as 



(G.7) 


(G.8) 


Using Equations (G.4 and G,5) the real part of the acoustic admittance. 


which is the quantity of interest as discussed in Section 4. 2. 5. 3, is 




L 


I J- + ™sl,I 

* -nm 

o 


(G.9) 


The real part of the acoustic admittance discussed in 
Section 4.2.5. 3 was determined from Equation (G.9) evaluated at n 
equal to infinity. The parameters in Equation (G.9) were taken from 
the numerical solution of the first order problem except for large 
ri where the asymptotic representations of the parameters were used. 



APPENDIX H 


EXPERIMENTAL DATA 
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Table 5 


Hydrazine Steady Strand Burning Rates 



Pressure 

(atm) 

Liquid Purity 
(%) 

Tube Inside Diameter 
(mm) 

2 

Burning Rate x 10 
(cm/sec) 


0.36 

99.4 

12.0 

1.61 

0.72 

99,4 

12.0 

2.01 

1.00 

99.4 

12.0 

2.42 

8.10 

99.4 

12.0 

8.74 

13.8 

99.4 

12.0 

16.1 

25.0 

99.4 

12.0 

26.2 

39.1 

99.4 

12.0 

44.5 

0.33 

98.6 

12.0 

1.55 

0.40 

98.6 

12.0 

1.60 

1.00 

98.6 

12.0 

2.40 

2.72 

98.6 

12,0 

4.25 

3.40 

98.6 

12.0 

5.00 

5.44 

98.6 

12.0 

6.50 

6.70 

98.6 

12.0 

8.61 

12.9 

98.6 

12.0 

16.72 

20.4 

98.6 

12.0 

25.8 

28.9 

98.6 

12.0 

35.5 

39.1 

98.6 

12.0 

48.2 

0.40 

98.6 

8.0 

1.72 

1.00 

98.6 

8.0 

2.55 

3.40 

98.6 

8.0 

5.15 

6.70 

98.6 

8.0 

8.81 

12.9 

98.6 

8.0 

17.3 

20.4 

98.6 

8.0 

25.9 

0.40 

98.6 

4.0 

2.20 

1.00 

98.6 

4.0 

3.25 

3.40 

98.6 

4.0 

5.95 

6.70 

98.6 

4.0 

9.80 

12.9 

98.6 

4.0 

17.5 

20.4 

98.6 

4.0 

26.0 

0.34 

95.6 

12.0 

1.41 

1.00 

95.6 

12.0 

1.89 

3.41 

95.6 

12.0 

3.65 

11.6 

95.6 

12.0 

6.30 

37.8 

95.6 

12.0 

12.5 

0.37 

92.6 

12.0 

1.25 

1.00 

92.6 

12.0 

1.70 

4.05 

92.6 

12.0 

3.30 

17.0 

92.6 

12.0 

6.10 

36.8 

92.6 

12.0 

10.6 
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Table 6 


Hydrazine Liquid Surface Temperatures 


Pressure 

(atm) 

Surface Temperature 
(K) 

0.51 

351.8 

0.82 

362.1 

1.00 

367.0 

3.37 

407.5 

6.60 

426.2 

10.1 

440.3 

13.0 

458.7 

17.2 

469.0 

19.8 

470.3 


Table 7 


Hydrazine Liquid Surface Oscillation Amplitude and Mean Burning 
Rate at a Mean Pressure of 1*54 atm and a Frequency of 1*32 Hz 

for Various Pressure Amplitudes 



£ 


V* X 10^ 

(cm/sec) 


0.051 

0.077 

0.311 

0.103 

0.075 

0.309 

0.153 

0.076 

0.311 

0.198 

0,074 

0.312 

0.245 

0.080 

0.315 

0.298 

0.086 

0.329 

0.335 

0.089 

0.335 
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Table 8 


Hydrazine Liquid Surface Oscillation Amplitude and Phase Angle at a 
Mean Pressure of 9.77 atm for Various Frequencies 


(0* 

(Hz) 

x|i/e 

(mm) 

$ 

(<“) 

0.52 

0.445 

98.1 

1.05 

0.201 

104.2 

2.88 

0.102 

107.4 

5.19 

0.065 

107.6 


Table 9 


Hydrazine Liquid Surface Oscillation Amplitude 
Mean Pressure of 1.18 atm for Various 


and Phase Angle at a 
Frequencies 


w* 

(Hz) 

x*i/e 

(mm) 

$ 

(<°) 

0.40 

0.121 

80.1 

1.01 

0.068 

59.8 

3.12 

0.017 

37.1 

4.85 

0.008 

-4.4 



